Emergence of a new pair-coherent phase in many-body quenches of repulsive bosons 
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We investigate the dynamical mode population statistics and associated first- and second-order 
coherence of an interacting bosonic two-mode model when the pair-exchange coupling is quenched 
from negative to positive values. It is shown that for moderately rapid second-order transitions, 
a new pair-coherent phase emerges on the positive coupling side in an excited state, which is not 
fragmented as the ground-state single-particle density matrix would prescribe it to be. 
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The fundamental question whether an interacting sys- 
tem of bosons is a single or fragmented condensate [lj- 
Q] has conventionally been answered in favor of the sin- 
gle Bosc-Einstein condensate. This conclusion is gained 
from highly symmetric (e.g. rotationally invariant) situ- 
ations, for example in antiferromagnetic spinor conden- 
sates H or rotating gases |f|. The fragmented conden- 
sate states occurring for these symmetrical Hamiltonians 
are extremely fragile against decay into a coherent state 
(a single condensate) by small perturbations coupling 
the single-particle modes Q, breaking, e.g., rotational 
symmetry. As a consequence, the fragmented states are 
also extremely sensitive to small time-dependent per- 
turbations like number fluctuations and generally exci- 
tations above the fragmented condensate ground state. 
Therefore, the instability of fragmentation towards a sin- 
gle condensate was traditionally assumed to be generic. 
However, a recent result within the two-mode approxima- 
tion for scalar bosons with general interaction couplings 
(that is when one moves away from points of high symme- 
try), has shown that a continous variety of ground state 
fragmentation can be obtained in a single trap 0. This 
type of fragmentation is, importantly, robust to pertur- 
bations coupling on the single-particle level. 

In the following, we investigate whether the robustness 
of continous fragmentation also persists against rapid 
changes on the dynamical many-body level, i.e. when in- 
teraction couplings rapidly change. Due to their many- 
body origin which is relying on the (relative) values of 
the couplings, fragmented states might then be expected 
to be more fragile. Here, we show that the sensitivity of 
fragmentation strongly depends on the range covered by 
a given coupling sweep. For a second-order dynamical 
quantum phase transition from the coherent to the frag- 
mented phase 0, the creation of single-trap fragmenta- 
tion from a coherent state is possible only for very slow 
sweeps of the coupling of bosonic pair exchange essen- 
tially determining the character of the ground state (sin- 
gle versus fragmented condensate [8(). For moderately 
rapid exponential sweeps, the final degree of fragmenta- 
tion is suppressed. The single-particle coherence measure 
gi (when averaged over time) vanishes after the sweep as 



well, while an analogously defined pair- coherence mea- 
sure 52 is macroscopically large. We therefore obtain a 
new pair-coherent phase which shows no single-particle 
fragmentation, in contrast to the ground state. The 
emergence of this phase is due to crossing the singu- 
lar point of vanishing pair-exchange coupling. On the 
other hand, for sweeps entirely on the positive exchange- 
coupling side, no excited-state suppression is obtained 
even for large sweep rates, and the degree of fragmenta- 
tion remains close to its ground-state value. We stress 
that these phenomena cannot occur in the fragmented 
phase of a conventional double- well with no pair exchange 
between sites included [Ioj |. 

Let us consider the following two-mode Hamiltonian 
describing the dynamics of particles in an arbitrary trap 
and for a two-body interaction Q 
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The Hamiltonian may be understood to result from go- 
ing one step beyond the familiar scmiclassical Gross- 
Pitaevskii theory i.e., by inserting the field opera- 
tor decomposition if! = ao'J'o + Si^i H — • into the full 
second-quantized Hamiltonian and truncating the expan- 
sion after the first two modes. We thus neglect 0(1) 
fluctuations on top of these (now generally still quan- 
tum) modes, which are assumed to be the only modes 
macroscopically populated with a finite fraction of the 
number of particles, TV; = O(N). The modes and two- 
body interactions are inhomogeneous and anisotropic in 
an independent way, leading to a set of interaction cou- 
plings {Ai} with essentially arbitrary relative magnitudes 
and signs. We assume that the modes have been de- 
termined, e .g., by the multiorbital mean-field method de- 
lineated in and integrated out. 

The sign of the pair-exchange coupling constant A3 
basically decides upon the classes of many-body ground 
states which can be obtained from the two-mode model 
(fTl) [8|1 . In contrast to optical lattices with pair-hopping 
[IJ-IlJ], here both the couplings A 3 and A± have gener- 
ally equal importance as A\ and A%. The Hamiltonian 
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(0 then describes, e.g., single-trap fragmentation in har- 
monically trapped quasi- ID and quasi- 2D gases in the 



weakly confining directions 15 1 



In what follows, we assume for simplicity that the 
single-particle states are degenerate, eo = ei, so that we 
can omit the first term in ([1]); this does not change our 
results in any essential way. We employ the following 
ansatz for the many-body wave function 
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which represents a linear superposition of Fock states 
\N — 1,1) = \l). We are interested in the time-dependence 
of the occupation amplitudes ipi(t), which comprise the 
information on the quantum many-body dynamics. At 
any instant, IV'iWI 2 gives the probability that I particles 
are in the second mode and N — I in the first. Both the 
amplitude and phase of the ?pi will influence the tem- 
poral behavior of the first- and second-order coherence 
functions g\ and gi, to be described below. 

The degree of fragmentation of an interacting many- 
body system, when maximally two field operator modes 
are macroscopically occupied, can be defined in an in- 
variant manner from the difference of the (macroscopic) 
eigenvalues of the single-particle density matrix [l], \ 
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Here, Ni — (hi) are the diagonal elements of the single- 
particle density matrix, while the first order coherence 
gi = ^(clqCii + a\ao) contains the real part of the off- 
diagonal correlator. For a coherent state with relative 
phase Ai9, g\ = \JN$N\ cos [A??], and J- vanishes. 

We assume that the interaction couplings Ai are large 
enough for the system to attain a fragmented state [15| . 
and that the couplings fulfill the condition A\ + A2 + 
2\A 3 \ - Ai > V t Q. To set the stage for the dy- 
namical case and introduce some notions, we first dis- 
cuss the stationary ground state solution for the ipi am- 
plitudes. Due to structure of the Hamiltonian, even 
and odd I sectors decouple from each other for A3 7^ 
(cf. Eq. (0 below), and one gets a solution of the ma- 
trix equations for ipi with a state for even and odd I 
separately, \<f>) = £ eveni 0;|Z), |$) - Eodd J *'!*>• Thcsc 
two states are degenerate in the large N limit. Conse- 
quently, any complex choice of a and b in the superposi- 
tion I*) = a \4>) + b |$) , \a\ 2 + \b\ 2 = 1 leads to a ground 
state at the same energy, with vanishing first-order coher- 
ence (71 and generally nonvanishing fragmentation T on 
the A3 > side. There results, for A3 > 0, a purely imag- 
inary (or vanishing) correlator (ajai) = ^2,u>iipfipi+i, 
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with u>i = yj (I + 1)(N — I) and we defined the relative 
phase between even and odd sectors 9 = arg(fe) — arg(a). 
The second line is valid in the continuum approximation 
of slowly varying occupation amplitudes, |$;| ~ |</>z±i|, 
and cui ~ W/±i. Hence gi, the real part of the above 
correlator, vanishes in the continuum approximation. 

The degree of fragmentation ([3]), as a function of \a\, 
\b\, and the relative phase 9 between a and b (which were 
assumed to be real in M) then evaluates to 
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We used the Gaussian ground state distribution of ~ 
(yJir/Aa)~ 1 / 2 exp[— (j - &)/(2a 2 )] representing 



m±i\ 

a solution in the continuum limit [4, 8]. We defined j = 
I - f , 6 = f - N Q = iVi - f as the shift from a 
maximally fragmented state and the distribution width 
fj - ^ 3 1 



iVEV4 with R{{A . }) = Ai+A2+2lM _ A4 
6 < N/2). For 6 = n/2 + 5 (\8\ < 1) and \a\ 



(when 
|6| 



1 / \/2, the degree of fragmentation is suppressed to T = 
>/R/N + <5 2 /2, while being maximal for a given set {Ai} 
at 9 = (mod tt). Adding a small (e.g. Cl = 0(A 1 /N)) 
Josephson type perturbation 
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on the A3 < side, the system is driven to a coherent 
state with a = b and T = 0. On the A3 > side, 
the ground-state fragmentation is insensitive to fi, the 
correction to T being of order Vl 2 /(N A3) 2 , provided the 
singular region around A3 = is excluded. 

We now proceed to the dynamical case, using A\ = 1 
as unit of energy in the following. We assume an expo- 
nential sweep of the form A3 (t) — (Ai— Af) exp[— at]+Af 
with A t = A 3 (t = 0) and Af = A 3 (t -> 00), and that all 
other couplings remain essentially constant during the 
sweep d t A 3 ^> <9 t Ai, <9 t A 2 , <9tAi [161]. The matrix equa- 
tions resulting from (0, (0 and (0 (h = 1) 
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are solved numerically, where the coefficients c/ = 
\A r (N - 1){N -l-l)+ \A 2 l(l - 1) + 7}A 4 (N - 1)1 and 
1 = ^(l + 2)(l + l)(N -Z-1)(JV -I). 

The Josephson parameter 17 is only dynamically ef- 
fective in a small window around A3 = according to 
the ratio £l/NA3(t) [cf. the ground state considerations 
above], and its effect is correspondingly small even for 
slow sweeps across the phase transition at A3 = 0; we 
illustrate this insensitivity by showing additional = 1 
data in Fig.[T] (green squares). While for second-order 
transitions, some small resonance peaks appear for slow 
sweeps [cf. Fig.[T](a)], at sufficiently large values of a, the 
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FIG. 1. Average asymptotic value of degree of fragmentation 
versus sweep rate a across A3 — 0. Red circles are from 
numerical solutions of (0 for SI = and solid lines a guide 
for the eye. For a second-order transition (a) A\ = Aa = 
1, A 2 = 0.5, At = A 3 (t = 0) = -0.2, Af = A 3 (oo) = 0.4, N = 
100, ground state final value of J-g S = 0.6, and for (weakly) 
first order in (b) they are A 2 = 0.6, TV = 200, J" gs = 0.71, 
AJ~ =1/3, others identical. The green squares are obtained 
by adding a perturbation of the form © with fi = 1. The 
inset shows J-(t) for a = 0.01; the blue arrow indicates where 
As = 0. Note the different vertical axis origin for (a) and (b). 



identical complete suppression of fragmentation occurs. 
We have verified that taking ft — 1 also does not quali- 
tatively alter the other = results to follow. 

During the sweep, the relative phase between the 
even and odd sectors, 9 = X^o{ ar S\ ! /'2fc+i) _ 
arg(V>2fc)}, V |-0i| 2 > 0> becomes time dependent. We will 
see that this has a crucial influence on the final degree 
of fragmentation T for sweeps from the coherent A3 < 
to the fragmented A3 > side. In the latter quantum 
phase transition, the singular point A3 = is crossed. 
At this point of vanishing pair-exchange coupling, the 
matrix problem in ([7|) becomes diagonal in the \l) Fock- 
space, and thus is easily solved. One Fock state \l) is ob- 
tained, with a fragmentation jump A J 7 = 1— | ^4^4r^ I ■ 
The transition is therefore generally of first order, with a 
discontinous change of fragmentation A J 7 at A3 = 0. We 
explore here both the second- and first-order cases of the 
dynamical quantum phase transition 17J when crossing 
the singular ^/-distribution point at A3 = 0. 



We observe a strong dependence on there being a 
ground-state fragmentation jump at A3 = 0, i.e. on 
the transition being second (AT 7 = 0) or first order 
(A J 7 7^ 0). The final average degree of fragmentation 
is defined as .F(oo) = t _ t J t ^ F{t')dt\ where ti ^s> t\ 
are two late times well after the phase transition point 
is crossed. For a 1/N, the final degree of fragmenta- 
tion .F(oo) quickly tends to zero in the second-order case 
(FiglH(a)). On the other hand, we obtain that JF(oo) de- 
cays much less rapidly with sweep rate a in the first-order 
case, cf. FigQ](b). We show the dependence of J"(oo) on 
AJ 7 in the left plot of Fig. [2] For the second-order case 
(a) in FigHJ the phase difference between the even and 
odd I sectors is rapidly driven towards an average value 
9 ~ 7r/2 after crossing A3 = for intermediate values 
of a (cf. Fig. [5]). The phase difference of 90 degrees bc- 



FIG. 2. Left: The final degree of fragmentation for a — 0.1 as 
a function of the jump at A3 = 0; N — 200 and the variation 
of A/7 7 is achieved by varying A2 from i to 1. Right: The 
average phase difference at late times versus AJ-". The inset 
shows that in the second-order case AF = 0, the average 9 
approaches n/2 at late times. Other parameters are identical 
with those of Fig.Q] £1 = 0. 



tween even and odd / sectors explains the strong suppres- 
sion of fragmentation in this regime of a 1, while for 
first-order transitions, 9 well after the transition is signifi- 
cantly less than n/2. The reason for this intimate relation 
between .F(oo) and 8(00) is that the argument leading to 
the ground-state result Q is still approximately valid 
when the condition arg(^ +1 ) — arg(^_i) ~ 7r(mod27r) 
is fulfilled for most I values, a property which we have 
verified. 

The above behavior for the phase-transition sweep 
needs to be contrasted with the case A 3 strictly positive 
during the whole sweep, imposing that A 3 3> 0(1/N). 
The average degree of fragmentation is not suppressed, 
even when a» 1 (and thus larger than the interaction 
couplings Ai). We illustrate this in Fig. [3] cf. the frag- 
mentation suppression observed for the quench case dis- 
played in Fig.Q] The insensitivity of fragmentation there- 
fore persists for dynamical changes of A3, provided the 
singular region around A 3 = is not traversed. 

We now discuss the pair-exchange coherence, measured 
by the expectation value 
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FIG. 3. Degree of fragmentation for a sweep entirely on the 
A3 > side for increasing sweep rate, where Ai = A4 = 
1,A2 = 0.5,7V = 200, Q, = 0. Even for large sweep rates of 
order Ai , J-(t) remains close to its ground state value, as well 
as 9 close its initial value 7r/4; \a\ = |6| = 1/V2 is constant. 
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a single-particle fragmented state, with the microscopic 
origin that the average phase difference between even and 
odd sectors of the many-body amplitudes approaches 7r/2 
after the quench. Rapid temporal changes of the interac- 
tion couplings of a many-body system can thus result in 
an emergent quantum phase with coherence properties 
strikingly different from those of the adiabatic ground 
state. 

This research was supported by the Brain Korea BK21 
program and the NRF of Korea, grant No. 2010-0013103. 



FIG. 4. Dynamical behavior of first- and second-order pair 
coherence gi(t) (top) and 32 (t) (bottom) for three different 
sweep rates a. In the first-order case (b), the pair coherence 
decreases for larger a while it remains on average identical 
in the second-order case (a). Parameter sets are identical to 
those of Fig.Q] with Q = 0. 



We observe, hrst of all, that the quantity 172 can be 
macroscopic, i.e., gi ~ 0{N 2 ), when the many-body 
state is neither coherent nor fragmented. In the ground 
state, it can be shown by direct calculation that for ei- 
ther sign of A3, there is pair coherence except right at 
A 3 = (where gi crosses zero), and that 32 is inde- 
pendent of the weights a and b. In particular, we have 
5 2 =i -N 2 /4 + n 2 /(8N 2 A 2 ) + a 2 /2 + 6 2 for A 3 > 0. 

The dynamical behavior of 172 across the sweep, along 
with that of its single-particle counterpart g\ is shown 
in Fig. 01 Whereas the degree of fragmentation is sup- 
pressed (see also Fig.QJ, and the time-averaged value 
of first-order coherence is zero, an essentially stationary 
pair condensate emerges, cf. the bottom row of Fig. 2] 
For a (weakly) first-order transition (in the case shown, 
A J 7 = 1/3), fluctuations become stronger and there is an 
increasing suppression of pair-coherence for larger a. 

While on the single-particle coherent side, g\ ^ 0, pos- 
itive pair-exchange coherence is trivially achieved due to 
the existence of first-order coherence, on the A 3 > 
side pair-coherence without single-particle fragmentation 
is the manifestation of a new pair-correlated phase emerg- 
ing after the sweep. We emphasize that the pair coher- 
ence does not result from attraction between bosons, as 
all A 
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are chosen positive on the A 3 > side 
We have shown that for moderatedly rapid variations 
of the pair-exchange coupling from negative to positive 
values in a second-order quantum phase transition, a new 
pair-coherent phase is created. In contrast to ground- 
state expectations, the resulting many-body state is not 
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